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Abstract 

In this paper, we argue that an extension of the Standard Model with a single 
leptoquark and three right-handed neutrinos can explain the excess in the first- 
generation leptoquark search at the LHC. We also find that when the leptoquark 
has similarly sized couplings to all three generations, it produces additional signals 
which will soon be tested in the second- and third-generation leptoquark searches, as 
well as in decay channels consisting of two mixed flavor leptons and two jets. If the 
leptoquark only couples to the first generation, on the other hand, two of the right- 
handed neutrinos need to be fairly degenerate in mass with the leptoquark while 
the other right-handed neutrinos mass should be much lighter. This hierarchical 
structure could explain dark matter and the baryon asymmetry of the Universe. 
These simple models may be regarded as benchmark models for explaining the 
excess, which can be tested in the next stage of the LHC running. 



1 Introduction 


Leptoquarks are particles that have both baryon number and lepton number, and appear 
in various beyond the standard model theories. For example, some grand unified theories 
(GUTs), such as SU(5) [1], SU(4)e <g) SU(2)i <g> SU(2)^ [2], or SO(IO) [3, 4] contain such 
particles. They are also found in technicolor models [5, 6] or other composite models 
[7]. Leptoquark searches, therefore, provide a powerful test for these theories beyond the 
Standard Model (SM). 

Recently, the CMS Collaboration announced that they observed an excess in the first- 
generation leptoquark search [8]. They analyzed the 8 TeV LHC data with an integrated 
luminosity of 19.6 fb -1 , and looked for first-generation leptoquarks in channels consisting 
of two electrons and at least two jets (eejj ), or an electron, a neutrino and at least two 
jets (eujj). An excess was found in both of these channels. However, this excess was only 
found after optimizing the event selection for leptoquarks with a mass of 650 GeV. The 
significance of the excess is estimated to be 2.4a and 2.6cr in the eejj and eujj channels 
respectively [8]. These excesses are, however, much smaller than those expected for a 
minimal extension of the SM with a leptoquark. In fact, the analysis carried out in Ref. [8] 
has excluded leptoquarks with masses less than 1005 (845) GeV if the branching fraction 
of the leptoquarks into an electron and a quark is equal to 1 (0.5). These constraints 
can be significantly relaxed, if the branching fraction to an electron and a jet is reduced. 
Actually, this suggests that leptoquarks have additional decay modes other than those 
with eejj and eujj final states. Not surprisingly, leptoquarks may have a richer structure 
than is found in these naive models. 

Stimulated by this observation, various interpretations of the excess have been pro¬ 
posed in the literature so far. For example, the authors in Ref. [9] introduce a coloron 
in addition to a leptoquark which could also explain the anomaly observed by the CMS 
Collaboration in the searches for a right-handed W gauge boson plus a heavy neutrino 
[10]. See also Refs. [11, 12] for an explanation of both of these excesses. Other leptoquark 
models with a dark matter (DM) candidate are considered in Refs. [13, 14], These mod¬ 
els may also explain the excess in the two opposite-sign same-flavor leptons with greater 
than two jets and missing transverse energy [15] channel. In supersymmetric theories, 
the existence of the A-parity violating interactions may induce the observed excess; this 
possibility is discussed in Refs. [16-18]. 

In this paper, we study another, in fact quite simple, model which can explain the 
excess in the first-generation leptoquark searches. This model consists of adding a single 
leptoquark and three right-handed neutrinos to the SM. In contrast to previous work, we 
only consider renormalizable interactions. As will be seen below, this model can actually 
account for the observed excess. It turns out that when the leptoquark couples to all 
three generations, signals in the second- and third-generation leptoquark searches are 
expected to appear in the near future. Furthermore, decay modes containing a pair of 
mixed flavor leptons and two jets are also expected. If the leptoquark only couples to the 
first-generation fermions, on the other hand, two of the right-handed neutrinos masses 
are required to be somewhat degenerate with the leptoquark mass. The remaining right- 
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Tabic 1: Quantum numbers of scalar lcptoquarks and possible renormalizable interactions 
containing leptoquarks. We follow the notation in Ref. [19]. 


Leptoquarks SU(3) C <g) SU(2) L <g) U(l) y Renormalizable couplings 


So 

( 3 , 1 , -1/3) 

So S 0 u^ R e^ R , SqQlQl, S 0 u R d R 

So 

( 3 , 1 , -4/3) 

So Squ r u r 

Si 

( 3 , 3 , -1/3) 

SiQ[4, S^lQl 

ct 

D 1/2 

( 3 , 2 ,+ 7 / 6 ) 

^\/-2^L e Ri ^i/2 U lt^ J L 

ct 

°l/2 

( 3 , 2 ,+ 1 / 6 ) 

Sjt/2 ^ r L l 


handed neutrino should be lighter than the electroweak scale. This hierarchical structure 
in the right-handed neutrino sector is quite interesting since the lightest right-handed 
neutrino can be a DM candidate while the other two could help explain the baryon 
asymmetry of the Universe (BAU). We will briefly discuss these possibilities. 

2 Model 

To begin with, let us describe our model with a 650 GeV leptoquark, which could account 
for the excess recently observed by the CMS Collaboration [8]. A comprehensive list of 
leptoquarks 1 is found in Refs. [20, 21]. In this paper, we focus on scalar leptoquarks since 
massive vector particles require additional complexity in order to explain their masses, 
i.e., the Higgs mechanism. We list scalar leptoquark models in Table 1. Among them, 
the first three types of leptoquarks cannot be assigned a definite baryon (B) and lepton 
(L) number such that the interactions of the leptoquarks conserve these charges [22], 
since they have both leptoquark and diquark couplings. In these cases, the tree-level 
exchange of leptoquarks induces proton decay [22-25] , which is stringently constrained by 
experiments . 2 For this reason, we do not consider these possibilities in this paper. For the 
other two leptoquarks, S [/ 2 and S^ 2 , there exist baryon-number violating dimension-five 
operators which are again problematic for proton decay. As discussed in Ref. [25], these 
operators can be forbidden by a Z 3 symmetry under which each held transforms as "0 —)• 
exp[27rf (B—L)/3]ip. This symmetry could be a remnant of the U(l) b-l symmetry [26-28], 
which can be realized naturally in GUTs . 3 Notice that this Z 3 symmetry cannot suppress 
the tree-level proton decay for So, So, and Si, since the renormalizable interactions given 

1 For a review, see the “LEPTOQUARKS” section written by S. Rolli and M. Tanabashi in Ref. [19]. 

2 In the case of So, additional W-boson exchange is required for proton decay to occur, since the 
SqUrUr interaction should include either a charm or top quark due to the antisymmetry in the color 
indices. The proton decay bound is still severe in this case. 

3 For instance, the vacuum expectation value (VEV) of a 672 of SO(IO) spontaneously breaks the 
U(1 )b_l symmetry into the Z 3 symmetry [29, 30]. 
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in Table 1 conserve B — L. In what follows, we will assume that the leptoquarks S\, 2 and 

S\ /2 do not lead to proton decay problems and focus on these two possibilities. 

As is often the case for a new scalar particle coupling to quarks and leptons, the 
addition of a leptoquark will, in general, cause flavor/CP problems. New interactions 
of a leptoquark with quarks and leptons may induce flavor-changing neutral currents 
(FCNC) and/or charged lepton flavor violation (CLFV), which are severely restricted by 
low-energy precision experiments. Furthermore, leptoquarks can also contribute to the 
electron and muon anomalous magnetic dipole moments, and the electric dipole moments 
of quarks and leptons; again the contribution should be small enough to be consistent with 
current experiments. For the existing constraints on leptoquarks from low-energy precision 
experiments, see Refs. [31-40]. These bounds are especially severe for leptoquarks which 
couple to both left- and right-handed quarks simultaneously, like S\ J , 2 in Table 1 [31]. For 

this reason, S^, is more promising. Although it may be possible to construct a viable 

model with S\j 2 , in this paper, we focus on S\, 2 , and for brevity denote it by S in the 
following discussion. 

As seen in Ref. [ 8 ], the CMS Collaboration observed slight excesses in both the eejj 
and the eujj channels. It may seem problematic that we consider the leptoquark 
since it can only give rise to the eejj or uujj modes and not the eujj mode. Further¬ 
more, the result in Ref. [8] indicates that leptoquarks have additional decay modes beyond 
those which produce the observed eejj and eujj final states. This motivates us to con¬ 
sider the possibility that the leptoquarks couple to other particles that are singlet under 
SU(3)c <8 U(1)em- With these additional couplings, it is possible to obtain the eujj -like 
events from decays to this singlet particle which yields a missing transverse energy sig¬ 
nal. In addition, these singlet particles allow us to hide the signature of leptoquarks if 
their masses are somewhat degenerate with the leptoquark mass, since the jets produced 
from the leptoquark decays become too soft to be detected. These additional hidden 
decay modes also result in a reduction of the branching fraction of the leptoquark to 
other states. Here, we will consider a singlet fermion which is coupled to the leptoquark 
through renormalizable interactions only . 4 The additional interactions of the leptoquarks 
are restricted by the SM gauge symmetries. The SM gauge symmetries indicate that these 
interactions should contain a quark held since leptoquarks are colored. The singlet fields 
should be fermionic in order to form a renormalizable Lorentz scalar from a leptoquark, 
a quark, and a singlet held. As we will see shortly, the leptoquark S actually has such an 
interaction with a singlet . 5 

The leptoquark S has B — 1/3 and L = — 1 in order for the model to preserve B 
and L. Thus, unless the leptoquark coupling with the singlet violates lepton number, 
the singlet fermions should have \L\ = 1. This charge assignment implies that these 

4 In Refs. [13, 14], on the other hand, leptoquarks are assumed to interact with a DM candidate through 
non-renormalizable interactions. 

5 Notice that ,S'j/ 2 does not couple to singlet fields via renormalizable interactions. In this sense, S\, 2 
is the only candidate for our considerations. However, it does still offer the possibility of hidden decays 
to other generations. 
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singlet fermions can be identified as right-handed neutrinos v R . Indeed the presence of 
right-handed neutrinos is well motivated by the observation of neutrino oscillations, which 
indicates neutrinos are massive. Based on the above discussion, we consider an extension 
of the SM to which a leptoquark S and three right-handed neutrinos v R . {% = 1,2,3) are 
added. 

The relevant interactions for our discussion are then given by: 


Ant = - A ~ K^Q^PrVjS^ 
- yipmMPrXLjY'H 3 - ^M,TPP R v t + h.c , 


( 1 ) 


where i,j are the generation indices; a represents the color index; ct, f3 are SU(2 )l indices; 
c indicates the charge conjugation; e a p is the antisymmetric tensor with 612 = —621 = + 1 ; 
Pl/r = (1 =F 75 )/ 2 . For later use, we write the SU( 2 )l components of the leptoquark S aa 
as 


S a = 



( 2 ) 


The relative size of the neutrino Yukawa couplings yij and Majorana masses are chosen 
such that the left-handed neutrinos have small masses consistent with the oscillation 
experiments. All of the interactions except for the Majorana mass terms conserve both 
baryon and lepton number, and thus do not induce rapid proton decay. The off-diagonal 
components of \j are strongly restricted by flavor experiments, and thus we assume the 
new couplings introduced above are diagonal, i.e., \j = A*5^- and = h t 5 ir Of course, 
this assumption can be relaxed within the experimental limits. 


3 Leptoquark signature at LHC 


At the LHC, leptoquarks are produced through gluon fusion and quark-antiquark anni¬ 
hilation. Their production cross sections are predominantly determined by the strong 
interactions, and the new coupling constants A* and A give only a negligible contribution 
to the leptoquark production since they are assumed to be very small. The reason for this 
will be clarified below. The NLO computation of the leptoquark pair production cross 
section is given in Ref. [41], 

After being produced, S u decays into the d R e c L or ulv r final state, while S d decays 
into the d R u R or d R v c R channel. The partial decay widths of these decay modes are as 
follows: 


nSu d m e c Li ) = T(S d -> d Ri u c Li ) 
T(S U ->■ u Li u Ri ) = T(S d ->■ d Li v Ri ) 


|Aj| 2 M L Q 

167T 


| hi 

2 m lq 

Ml 1 

1 -tli 

167T 

m Iq . 


(3) 


where Mlq and M Ri are the masses of the leptoquark S and right-handed neutrinos, 
respectively. A pair production of S u can yield both the eejj and the evjj signals, while 
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Sd gives uujj signals only. Since the first-generation leptoquark searches by the CMS 
Collaboration only look for the former two signals, the production of Sd is irrelevant for 
our discussion. 

In these searches, the event selection for both eejj and uejj signals is optimized for 
each leptoquark mass by imposing additional cuts beyond those referred to as the prese¬ 
lection cuts. A list of these cuts can be found in Ref. [8]. After all the cuts are imposed, 
the number of the leptoquark events is estimated to be 125.85(58) and 37.22(37) in the 
eejj and eujj channels, respectively, for Mlq = 650 GeV. Since the NLO production cross 
section for a 650 GeV leptoquark is determined to be 13.2 fb [41], the signal acceptances 
for the leptoquarks are 48.6% and 28.8% for eejj and eujj [9] channels, respectively, 
for the integrated luminosity 19.6 fb -1 . On the other hand, 36 events are observed with 
20.49 ± 2.14(stat) ± 2.45(syst) events expected in the eejj, and in the eujj channel 18 
events are detected with 7.54 ± 1.20(stat) ± 1.07(syst) events expected. Therefore, to 
explain these excesses the branching fractions of the S u —> de c and S u —* uu c are required 
to be 


BR(S' U —> de c ) ~ 0.35 , 

BR (S u -A uu c ) ~ 0.20 . (4) 


Clearly, the leptoquarks must have additional decay channels. 

One of the most straightforward and simple explanations of the additional decay chan¬ 
nels is that the leptoquark also decays into second- and third-generations quarks and 
leptons. We refer to this scenario as Model I in what follows . 6 For the assumption 
Mji r <C Mlq, the branching fractions of S u are given by 


BR (S u ->• de c ) 
BR (S u ->• uu c ) 


|Ai | 2 

Yli |Ai| 2 + Yli \hi\ 2 

ElM 

| A *| 2 + l ^*| 2 


( 5 ) 


Therefore, if we take the appropriate size for the couplings A 2 and A 3 , we are able to explain 
the preferred branching fractions found in Eq. (4). The decay modes of the leptoquarks 
into the second- and third-generation quarks and leptons are restricted by searches for 
leptoquarks which couple to only the second- or third-generation quarks and leptons. For 
second-generation leptoquarks, the CMS Collaboration studied the final states consisting 
of two muons and at least two jets (the /i/ijj channel) and one muon, and at least two jets, 
and missing transverse energy (the fiujj channel) with the 8 TeV 19.6 fb -1 data [42], The 
former channel leads to the constraint BR(S' U —y dfi c ) < 0.27 for a leptoquark of mass 650 
GeV and the latter channel leads to BR(S' U —> d/A) < 0.64 assuming BR(S' U —» uu c ) = 0.20 
for the same leptoquark mass. They also searched for a third-generation leptoquark with 

6 As already mentioned in footnote 5, Sy 2 given in Table 1 may also explain the CMS excess in a 
similar manner, though this type of leptoquark in general suffers from severe low-energy experimental 
constraints. 
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(a) BR(S'„ —>• fis) vs. BR(S7 —► rb) 


(b) IV jj channels 


Figure 1: (a) Branching fractions to the second and third generations that can explain the 
excess in Ref [8] are plotted in a red solid line with the corresponding band showing the 
uncertainty coming from the experimental errors in the event rates. Right gray-shaded 
region is excluded by the second-generation leptoquark search in the /i/ijj channel [42], 
while the upper blue shaded region is disfavored by the third-generation leptoquark search 
[43]. (b) Scattering cross sections for the IV jj modes as functions of BR(S' U —» /is). Red 
solid, blue dashed, and green dotted lines show the cross sections of the e/ijj, erbj, and 
firbj modes, respectively, with the corresponding bands indicating the uncertainties from 
the error in the observed event rates. The gray shaded region in the right side is excluded 
by the second-generation leptoquark search in the /j/ijJ channel [42], The leptoquark 
mass is set to be 650 GeV. 


the 8 TeV 19.7 fb -1 data set, and found the bound on the branching fraction for a 
leptoquark decaying to a tau lepton and a bottom quark to be BR(S' n —» rb) < 0.74 
[43]. The decay mode of a leptoquark to a top quark and a neutrino can in principle be 
constrained by the stop search, though the current constraint is still very weak [44, 45]. 
The ATLAS Collaboration also searched for leptoquarks of these three generations with 
their 7 TeV data set [46-48] , though currently their constraints are weaker than the CMS 
results. In Fig. 1(a), we show the preferred branching fractions to the second and third 
generations in a red solid line with the corresponding band representing the uncertainty 
coming from the experimental errors in event rates. Here, the right gray-shaded region 
is excluded by the second-generation leptoquark search in the /i/ijj channel [42], while 
the upper blue shaded region is disfavored by the third-generation leptoquark search 
[43]. We find that to explain the excess in Ref. [8] the leptoquark must have a sizable 
branching fraction into the third-generation fermions, BR(S' U —> rb) > 0.2. This rather 
large branching fraction may be probed in the upgraded LHC experiments. Note that 
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since only the branching ratios are relevant to our discussion, the size of the new couplings 
can be small. Thus, we need not worry about non-perturbativity of the couplings nor the 
constraints coming from single leptoquark production at HERA experiments [49, 50] and 
the LHC experiments [51-53]. Because the leptoquarks can decay to second- and third- 
generation quarks and leptons, this model can be verified by searches for leptoquarks 
with final states consisting of different flavors of charged leptons and two jets (the IV jj 
channels with l ^ V, and V,l = e,/i,r). In Fig. 1(b), we plot the scattering cross sections 
for these final states as functions of BR(S' n —> /is). Here, the red solid, blue dashed, and 
green dotted lines show the cross sections of the e/ijj, erbj, and firbj modes, respectively, 
with the corresponding bands indicating the uncertainties from the error in the observed 
event rates in Ref. [8] . Again, the gray shaded region on the right side is excluded by the 
second-generation leptoquark searches in the nnjj channel [42], The leptoquark mass is 
set to be 650 GeV. This figure shows that the cross sections for the IVjj channels can be 
several fb, which are within the reach of the LHC experiments. 

Here, we note that although we have assumed M Ri -C M L q in this analysis, some or 
all of the right-handed neutrinos can have masses of (9(100) GeV without modifying the 
results too much. This freedom is important when we discuss the possibility of right- 
handed neutrino DM, as will be seen in Sec. 4. 

Next, we consider the case where the leptoquark only couples to the first-generation 
fermions, i.e., X 2 = A 3 = 0. We call this scenario Model II. For this case, to realize the 
branching fractions in Eq. (4), we need to hide some of the decays into the right-handed 
neutrinos. To accomplish this, we assume that some of the right-handed neutrinos are 
degenerate in mass with the leptoquark. This makes the decays of the leptoquark to the 
degenerate right-handed neutrinos invisible since jets in this decay process become too soft 
to be detected. For example, the first-generation leptoquark searches in Ref. [ 8 ] require 
that the leading jet should have a transverse momentum larger than 125 GeV. Thus if 
the mass differences between the leptoquark and right-handed neutrinos are less than 
~ 100 GeV, the jet from the decay of the leptoquark evades the selection. From Eq. (4), 
we see that the invisible decays to right-handed neutrinos should have a branching ratio of 
about 0.45, which is twice as large as the visible neutrino branching ratio BR(S' n —> uu c ). 
This observation is in line with two of the three right-handed neutrinos being somewhat 
degenerate in mass with the leptoquark. In this case, each branching ratio is given by 

BR(S„ -> de‘) = -At- , 

N 2 + E,M 2 [i-b £] 2 


BR (S u -> uv c ) 


_ \hi \ 2 

I'M 2 + Yli IM 2 [1 


Mr. 

M, 


'ky 

LQ J 


) 


BR(invisible) = ^ 


I h s 


M| 


Ri 


j=2,3 |Al| 2 + EJ^| 2 [l- fe] 2 L 

L/CJ 


m lq. 


where we have assumed M Rl <C M R2 ,M Ra . This gives \hi\ ~ 0.76|Ai|. 


( 6 ) 


Moreover, if 
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|/i 2 | = \h 3 \ = h and M L q — M Ri = M L q — M Rs = 100 GeV, then h ~ 2.8|Ai|. At present, 
the coupling Ai is a free parameter. However, this prediction for the relation among 
the couplings can be important when the lightest right-handed neutrino is the dominant 
component of DM in the Universe, as we will discuss in the next section. Also, Ai will be 
constrained in this case. 

Finally, we consider constraints on these models from the clijet plus missing transverse 
energy (jj + Ifir) searches. After being produced, 100% (4%) of Sd ( S u ) pairs decay into 
two quarks and two neutrinos giving the jj + fix events in the detector. For a 650 GeV 
leptoquark decaying to jj + fir, the current upper limit on the cross section is about 
a few xlO fb [54, 55],' which is larger than the leptoquark production cross section. If 
some of the right-handed neutrinos have similar masses to that of the leptoquark, then the 
constraints are significantly relaxed. Anyways, we expect that the LHC run II experiments 
can observe events in this mode. This is a powerful consistency check of our models. 


4 Right-handed neutrino DM and the baryon asym¬ 
metry of the Universe 


In the previous section, it was found that the CMS measurement can be explained by a 
leptoquark by adjusting the branching fractions (couplings) of the leptoquark. However, 
the upper bound on the couplings is independent of the CMS signal and a rather weak up¬ 
per bound on the couplings is placed by HERA [49, 50] and CMS [51-53] single-produced 
leptoquark searches. In fact, order-one couplings are still allowed. The lower bound on 
the couplings of the leptoquark come from requiring that the leptoquark decay be prompt. 
This amounts to having a decay length shorter than about 1 mm. If the decay length is 
longer than 1 mm, then the constraint from searches for displaced vertices can be severe 
[56]. The total decay width for the leptoquark discussed above is 


r,s = = b.svi - y. 


I Aj| 2 + \hj\ 


1 - 


m Iq 


Af, 


LQ 


|Ai 


M 


LQ 


167T BR(S' n —y de c ) 167T 


(7) 


for both Models I and II. Then, we find the decay length to be 

'7.3 x 10" 8X 2 


cT^ 1 ~ 1 mm x 


650 GeV 


Mi 


LQ 


X 


Ai 


( 8 ) 


As can be seen from this expression, if the couplings are smaller than 10“' their decays 
can no longer be considered prompt. Finally, we find that the allowed range for the 
leptoquark couplings is 10 -7 < Ai < 1, and they are basically unrestricted by collider 
searches. 7 8 

7 The ATLAS limit given in Ref. [55] only shows the results obtained for a combined 2-6 jets plus 
missing energy search. This makes it rather difficult to ascertain the constraints on the production cross 
section for this mode alone. 

8 Low-energy precision measurements, on the other hand, give fairly strong constraints on the lepto¬ 
quark couplings. For example, measurements of the parity-violating transition in cesium atoms [57, 58] 










The size of the leptoquark couplings A* and hi, however, becomes important once we 
consider the phenomenology of right-handed neutrinos. Right-handed neutrinos may play 
a significant role in both particle physics and cosmology since they explain not only small 
neutrino masses via the seesaw mechanism [61], but also DM and the baryon asymmetry 
of the Universe. Indeed, it has been proposed that an extension of the SM with three 
right-handed neutrinos can account for all of these phenomena. This scenario [62] is 
called the v Minimal Standard Model (z/MSM). * * * * * * 9 It is therefore quite interesting to see if 
adding a leptoquark still gives a model which affords a DM candidate and can generate 
the desired baryon asymmetry. We will see that the range for the leptoquark couplings 
will be drastically reduced by the above requirements. 

4.1 Vanishing Leptoquark Couplings 

As a warm up, let us first consider the limit of vanishing hi and Aj. In this limit, the 
model considered in this article is exactly that of the z/MSM [62], for which DM and 
baryogenesis can be explained. Later, we will discuss the effects the leptoquark couplings 
have on the phenomenology of these models. As we will see below, these couplings can 
drastically alter the situation. In the z/MSM, the lightest right-handed neutrino is assumed 
to have a mass of 0(10) keV and to be the DM of the Universe. The masses of the other 
two heavy right-handed neutrinos are quasi-degenerate and are of 0(0.1-100) GeV in 
order to generate the baryon asymmetry. Both DM production and baryogenesis are 
non-equilibrium processes. For this to be the case, it is assumed that no right-handed 
neutrinos are generated during reheating. After reheating, the right-handed neutrinos 
are generated by various scatterings of the thermal bath. Since the oscillations among 
different generations of right-handed neutrinos can be a non-equilibrium process which 
violates lepton number, a lepton asymmetry is generated in the right-handed neutrinos 
among the different generations. The lepton sector includes extra CP phases and thus 
offers new sources for CP-violation, which is required to produce a baryon asymmetry. The 
asymmetry in the right-handed neutrinos is converted to an asymmetry in the left-handed 
neutrinos via oscillations. This asymmetry is then converted to a baryon asymmetry via 
the sphaleron [64], Importantly, the right-handed neutrinos do not reach equilibrium 
until after electroweak symmetry breaking, when the sphaleron shuts off, preventing the 
generated asymmetry in the right-handed neutrinos, and thus in the left-handed neutrinos, 
from being washed out. Because the sphaleron process is inert when the right-handed 
neutrinos reach equilibrium, the baryon asymmetry is preserved. 

The production mode for the right-handed neutrinos is predominantly through the 
scattering tt —>■ h* —> for the z/MSM, where N, t and h denote right-handed neutrinos, 

restrict the first-generation leptoquark couplings so that |Ai| < 0.22 [52] for a 650 GeV leptoquark. The 

Qweak experiment [59] , which determines the weak charge of the proton by measuring the parity-violating 

asymmetry in elastic electron-proton scattering, gives a similar constraint on this coupling. The most 

severe constraint for Model I comes from AT —> fi~e + decays. In fact, the constraint on the leptoquark 

couplings is |A 2 | < 0.9 x 10 -5 [52, 60]. As we will see, however, these are not the strongest constraints 

on the couplings of our model. 

9 For a review, see Ref. [63]. 
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top quark and the Higgs boson, respectively. When the temperature T is much higher 
than the masses of the particles participating in the scattering process, this scattering 
leads to a production rate of 

Tat = n EQ (a t i^ U N v) ~ l(T 3 |?/i i | 2 T , (9) 


where i ^ 1, n E q indicates the number density of a massless particle and (a t t^vNv) is 
the thermally averaged scattering cross section times relative velocities. In the i/MSM, 
the Yukawa couplings of the neutrinos with the Higgs are quite small in order to explain 
small neutrino masses. For the heavier neutrinos they are of order y l3 < 1CT 7 when i ^ 1 
while i/ij < 10~ n for i = l. 10 The heavier right-handed neutrinos will then come into 
equilibrium much earlier than the lightest right-handed neutrino. They will come into 
equilibrium when Tat ~ H (H is the Hubble expansion rate) which gives 

Teq “ ~ (li ^) 2 x o(10) GeV • (10) 

where Teq is the temperature at which the heavier neutrinos come into thermal equilib¬ 
rium; Mp = 1.22 x 10 19 GeV is the Planck mass; g* is the effective degrees of freedom 
for relativistic particles in the thermal bath. A more detailed calculation shows that to 
ensure right-handed neutrinos are out of thermal equilibrium before the sphaleron process 
shuts off, the Yukawa couplings should be smaller than 2 x 10~' [62], and thus the heav¬ 
ier right-handed neutrino masses should be Mp 2 3 < 20 GeV. Furthermore, this scenario 
means that the largest asymmetry in the right-handed neutrinos, and likewise in the left- 
handed neutrinos, is produced just before the sphaleron shuts off. This larger asymmetry 
can then offset the small mixing of the right- and left-handed neutrinos and produce the 
needed baryon asymmetry. 

DM in the ffiVISM is the lightest right-handed neutrino. It is well known that a sterile 
neutrino with a mass of 0(10) keV is a good warm DM candidate [65, 66]. The Yukawa 
couplings of the lightest right-handed neutrino, if it is DM, need to be sufficiently small 
so that its lifetime is longer than the age of the Universe. The dominant decay mode of 
the right-handed neutrino is the N —> ?>v channel. By evaluating the decay rate, one finds 
that an 0(10) keV right-handed neutrino naturally realizes a lifetime much longer than 
the age of the Universe. As it turns out, however, this is not the most severe restriction 
on the couplings of the lightest right-handed neutrino if it is DM. Its couplings are more 
severely constrained by searches for diffuse X/y rays [67-70]. If the decay width is too 
large, it will produce diffuse X/y rays which can be detected. To avoid these constraints, 
the lifetime needs to be longer than about 10 28 s which gives 


T/Vp—>71/ 
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(ii) 


10 In the presence of the leptoquark S, the Yukawa coupling can also be induced at the loop level 
since the leptoquark couples to both left- and right-handed neutrinos. This contribution is potentially 
dangerous since it could generate Yukawa couplings which are too large. However, it turns out that 
this contribution is smaller than that from the seesaw mechanism for the scenarios discussed below, i.e., 
leptoquark Yukawa couplings 0(1O -4 ). 
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where a is the fine-structure constant; Gf is the Fermi constant; v ~ 246 GeV is the 
Higgs VEV; N± denotes the lightest right-handed neutrino DM. As can be seen from these 
expressions, the decay is facilitated by a mixing of the left- and right-handed neutrinos 
and a lU-boson loop. By inhibiting the oscillations between the right- and left-handed 
neutrinos, this decay mode can be drastically suppressed. Suppressing oscillation amounts 
to taking smaller values of i/ij. 

Because the Yukawa couplings for the lightest right-handed neutrino are extremely 
small, its production from top scatterings is suppressed and therefore not enough DM is 
produced. This problem can be alleviated if a large asymmetry is generated in the left- 
handed neutrinos which persists to scales below electroweak symmetry breaking [71]. This 
large asymmetry in the left-handed neutrinos can be produced via neutrino oscillation of 
the two heavier right-handed neutrinos as they freeze out or from resonantly enhanced 
decays [72], For the oscillation production mechanism, the asymmetry is produced via 
mixing among the right-handed neutrinos. When the mass scale associated with mixing 
of the right-handed neutrinos, i.e., the differences of the right-handed neutrino masses, 
are similar in size to the Hubble parameter, an asymmetry develops among the different 
species of the right-handed neutrinos. This asymmetry is then converted to an asymmetry 
in the left-handed neutrinos via mixing between the right- and left-handed neutrinos. 
Since this oscillation occurs after freeze-out of the right-handed neutrinos, the left-handed 
neutrinos do not oscillate back to the heavier neutrinos. A similar mechanism is used to 
produce the baryon asymmetry as the right-handed neutrinos are being generated by 
scatterings of the thermal bath before equilibrium has been reached [62, 73]. However, 
the production of the baryon asymmetry must occur at a temperature above (9(100) 
GeV before the sphaleron shuts off, whereas the asymmetry needed for DM production 
is produced below the scale where the sphaleron shuts off. If the asymmetry in the left- 
handed neutrinos is produced via decays, the masses of the right-handed neutrinos still 
need to be degenerate in order to enhance the CP-violating parameter and produce a large 
enough asymmetry in the left-handed neutrinos [74], Once this asymmetry in the left- 
handed neutrinos is produced, the thermal effects associated with this non-zero chemical 
potential alter the neutrino mass matrix and so alter the mixing relations between the 
left- and right-handed neutrinos, like the Mikheyev-Smirnov-Wolfenstein effect [75]. In 
fact, for some temperature the mixing becomes of order one. If the temperature of the 
Universe for which the mixing is order one is still greater than the mass of the lightest 
right-handed neutrino, the left-handed neutrinos can be converted into the lightest right- 
handed neutrino giving a sufficient density to account for DM. 


4.2 Non-Vanishing Leptoquark Couplings 

For non-zero values of hij and i/ij, things change drastically since they couple the neutrinos 
with SM fields. These additional couplings provide additional scattering processes which 
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can generate the right-handed neutrinos. The most important process among them is the 
leptoquark-gluon scattering, gS —> QlN, which gives a production rate of order 


Ttv = n-EQ(cTgs^Q L Nv) ~ 10 2 a s \hij\ 2 T . 


(13) 


The equilibrium temperature for the right-handed neutrinos is then 


Teq ~ 10 a s | h 


Mr 


13 1.66 gl' 2 


10 x 
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(14) 


If the couplings are taken as small as is allowed by the prompt decay limit, we could 
get a situation similar to the z/MSM for baryogenesis. Although this may work, it does 
significantly restrict the parameter space of the model. 

However, this production mechanism is very problematic for DM in Model II. Since the 
LHC measurements require the coupling of each right-handed neutrino to be of similar 
size, the amount of the lightest right-handed neutrino produced will be similar to the 
other right-handed neutrinos. This corresponds to a near thermal abundance for the 
lightest right-handed neutrino, which will, in general, be too much DM . 11 Since its decay 
modes become severely inhibited for energies below the leptoquark mass, the right-handed 
neutrinos produced will stick around and overdose the Universe. 

Below, we will discuss how the lightest right-handed neutrino could still be DM in 
this scenario. However, first we discuss constraints on the couplings of the lightest right- 
handed neutrino to leptoquarks if it is indeed DM. These additional couplings will give 
additional contributions to T^-h^. Because the interaction of the leptoquark with the 
neutrinos are vector like, this decay has a contribution coming from a loop involving a 
leptoquark and a quark [76]. A lifetime of the lightest right-handed neutrino of order 
10 28 s, again required by the constraint coming from diffuse X/y rays [69, 70], requires 
parameters of order [76] 12 


10 _4 \ 2 /10 “ 4 \ 2 / M lq \ 5 / 4 keV 
A, ) V ^ J \650 GeV ) V M Rl 


x 10 28 s . 


(15) 


Fortunately, these couplings are much larger than those needed for the prompt decay of 
the right-handed neutrinos. If the couplings lie right at the upper edge allowed by this 
constraint, the leptoquarks with right-handed neutrino DM could explain the 3.5 keV line 
observed in X-ray spectra [77]. 

As can be seen from the previous paragraph, the constraints on the couplings of right- 
handed neutrino DM to the leptoquark are not so severe. The problems for producing 
the right-handed neutrino DM arise from thermal processes for T > Mlq. Below this 
temperature the leptoquarks decouple from the thermal bath and have little effect on 

11 In Model I, it may be possible to suppress production of the lightest right-handed neutrino by taking 

h\ <C /l2,3- 

12 This constraint on the couplings is more mild than that found in Ref. [76]. This is due to the lightest 
right-handed neutrino only coupling to the down-type quarks. 
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the cosmology. One way to circumvent this problem is, therefore, to have a very low 
reheat temperature. If the reheat temperature Tr is smaller than the leptoquark mass, 
the production of the lightest right-handed neutrinos from Sg NQ will be Boltzmann 
suppressed. This allows us to take larger values of A and h while keeping right-handed 
neutrinos from being over produced by the thermal bath. For larger values of h and A, 
the scattering process Ld —> QN will become important. Through this process, the right- 
handed neutrino DM can be non-thermally produced. This is the same as the non-thermal 
equilibrium DM scenario discussed in Refs. [30, 78]. In this scenario, the right-handed 
neutrino DM is only weakly coupled to the thermal bath. This is due to the small 
couplings A and h and because the mediator of this process, the leptoquark, has a mass 
somewhat larger than the reheating temperature. The lightest right-handed neutrinos 
are then produced by scatterings of the thermal bath. Because of the small couplings 
and large leptoquark mass, the right-handed neutrinos never reach thermal equilibrium. 
However, for particular values of the couplings, leptoquark mass, and reheat temperature 
enough of the lightest right-handed neutrinos can be produced to account for DM. 

If the right-handed neutrinos come into equilibrium, they will overdose the Universe 
and could cause problems for galaxy formation [65] . The right-handed neutrinos will come 
into equilibrium when the production rate is equal to the Hubble parameter. To guarantee 
that the lightest right-handed neutrino never comes into equilibrium, we calculate the 
production cross sections for T < T R < M L q to leading order 

3 'j' 2 

(VdL^QN, V) = V ) ~ > ( 16 ) 

i i L Q 


where Nq = 3 is the SU(3)c color factor and N w = 2 is for the doublet factor of SU( 2 )^. 
The right-handed neutrino DM is not in thermal equilibrium if the production rate is 
much smaller than the Hubble expansion rate. This condition, (<7Q7Vi^diT) n !v Q < H, is 
always satisfied if 


J] \\ihi\ 2 < 1.5 x 1(T 14 



( \ 
\650 GeVy ’ 


(17) 


where we have used the inequality T < Tr < Mlq. As can be seen, if A* ~ A ~ A for 
A i, hi 7 ^ 0, then A < 3.5 x 1CT 4 satisfies this condition. For this case, the right-handed 
neutrino DM cannot come into equilibrium after reheating. 

After reheating, the number density for the right-handed neutrinos is effectively zero. 
Since the number density is effectively zero, the annihilation rate of the right-handed 
neutrinos is negligible. The Boltzmann’s equation for generating DM can be simplified to 
include only the production part [30, 78] 


dY N (x ) 
dx 



MpMn 1 


i^dL^QN! v ) ^2 

X 2 EQ ’ 


(18) 


where x = M Rl /T, Yn = un/s, and Yeq = ueq/s; tin is the number density of the 
lightest right-handed neutrino DM; s is the entropy of the Universe; g s and g p are the 


13 







effective degrees of freedom for the entropy and energy density respectively; {cr d i^Q Nl v) is 
as above. By integrating the above equation, we obtain the required reheat temperature, 


Tn ~ 400 


/ DdM^ 2 \ 3 

V 0.12 ) 



( Mlq \ 
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/10 keV\ 

V M Rl J 




(19) 

where we set g s = g p = g*. Nq is the number of generations participating in the pro¬ 
duction process, which is model dependent. Qdm ^ 2 is the present DM density parameter. 
This reheat temperature falls exactly where is needed for Ni to be a good DM candidate. 
Namely, it is just below the lcptoquark mass so that the leptoquark interactions are sup¬ 
pressed. Furthermore, it is above the sphaleron process which means the non-equilibrium 
dynamics could generate a baryon asymmetry via the lepton asymmetry induced by os¬ 
cillations in the heavier right-handed neutrinos, as discussed above. This scenario is more 
plausible for Model I, where we may take the heavier right-handed neutrino masses to be 
0(1-10) GeV. The baryon asymmetry may also be generated through decays after the 
right-handed neutrinos freeze out. Again, the asymmetry in the left-handed neutrinos is 
converted to the baryon asymmetry through the sphaleron process. Since the constraints 
from diffuse X /7 rays only pertain to the lightest right-handed neutrino, the leptoquark 
couplings could also be adjusted to increase or decrease the baryon asymmetry as needed. 

To determine the precise relic density of the right-handed neutrino DM and the baryon 
asymmetry in our models, it is necessary to numerically solve a set of coupled Boltzmann 
equations with thermal effects adequately taken into account. This calculation is necessary 
to determine the leptoquark couplings with better accuracy. In particular, a more detailed 
computation is needed to determine if these models can explain the 3.5 keV X-ray line 
excess [77]. This analysis will be left for future work. 


5 Gauge coupling unification 


It is widely known that the SM is not compatible with the minimal SU(5) GUT [1], First, 
the SM gauge couplings approach each other but never unify sufficiently [79]. Second, 
the GUT scale is predicted to be too low to evade the current proton decay constraints. 
Therefore, the simplest SU(5) GUT is ruled out. 

Since the leptoquark S is charged under the SU(3)c < 8 > SU(2 )l <g) U(l)y gauge inter¬ 
actions, it can affect the running of the SM gauge couplings. Thus, the presence of the 
leptoquark may solve the problems mentioned above and revive the minimal SU(5). In¬ 
deed, the leptoquark we consider can be embedded in a 10, 15, 40, ... of SU(5). This 
allows us to construct a GUT model where S' is a component of an irreducible represen¬ 
tation of SU(5). 

At the two-loop level, the running of the SM gauge coupling constants g a (a = 1, 2, 3) 
is given by 


dg a 

^ djjL 


-J—b (1) g 3 + —^— 
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where (]\ = ^5/3g' and /*, (k = u , d , e) denotes the Yukawa matrices of the SM fermions. 
The coefficients b£\ , and c a fc in the SM are given in Ref. [80]. The change in the 

coefficient of the beta function from S is found to be 
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Here, we neglect the effects of the leptoquark Yukawa couplings since they are assumed 
to be very small. 

In Fig. 2(a), we plot the running of the SM gauge couplings with a leptoquark S 
in solid lines. The blue, green, and red lines show a~ l (a = 1,2,3), respectively, with 
a a = (y^/(47r). The dashed lines show the running due to SM fields alone. In this 
computation, we use the two-loop renormalization group equations given above. However, 
we have not included threshold corrections at the GUT scale since they are dependent 
on an unknown mass spectrum. From this figure, we see that although the presence of 
a leptoquark improves gauge coupling unification, the deviation is still sizable and the 
unification scale is still rather low. 

The situation is drastically improved, however, if there is an additional leptoquark 
around the TeV scale. A leptoquark with a mass of 0(1) TeV can easily evade the 
current LHC bounds. In Fig. 2(b), we show the gauge coupling running for this case. We 
also show the case where two leptoquarks and a Higgs doublet 13 are added to the SM in 
Fig. 2(c). This situation was already noticed in Ref. [81]. As can be clearly seen, gauge 
coupling unification in these two cases is far more precise. Accordingly, an extension of 
the SM with leptoquarks provides an alternative and promising framework for an SU(5) 
GUT without supersymmetry. 

In the above cases, the unification scale is predicted to be ~ 7 x 10 14 GeV, which is 
still too low if one considers the bound from proton decay. In the non-supersymmetric 
GUTs, proton decay is induced by the exchange of the heavy SU(5) gauge bosons. The 
dominant decay channel is then the p —> e + 7T° mode. The present experimental limit on 
this channel is r(p —> e + 7r°) > 1.4 x 10 34 years [82]. To evade this constraint, we need to 
assume that the heavy gauge boson mass M x is > 3 x 10 15 GeV, which is much heavier 
than the unification scale found above. With this disparity between scales, we expect 
sizable threshold corrections at the GUT scale and thus the gauge couplings do not need 
to precisely unify [83]. In this sense, all of the models in Fig. 2 can be compatible with the 
SU(5) GUT framework and thus they are quite promising. Anyways, all of these scenarios 

13 The contribution of a Higgs doublet to the beta-function coefficients is given by 
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(c) Two leptoquarks and one additional Higgs 

Figure 2: The running of the gauge couplings with (a) one leptoquark, (b) two leptoquarks 
or (c) two leptoquarks plus one additional Higgs doublet added to the SM are shown in 
solid lines. The SM running is shown in dashed lines for comparison. 


give a relatively light SU(5) gauge boson, and thus future proton decay experiments should 
detect proton decay in the p —y e + 7r° mode. This is a distinct prediction of these models 
compared with the SUSY SU(5) GUT scenario. In the SUSY GUTs, the dominant decay 
channel is generically the p —> K + u channel [84]. Additionally, since proton decay in 
our scenario occurs via the gauge interactions, the various decay modes are related by 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. Thus, this model will give specific 
predictions for ratios such as BR(p —» e + /l°)/BR(p —> e + 7r°). As a consequence of this, 
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the leptoquark extension of the SM discussed in this paper offers a simple way to achieve 
an SU(5) GUT model [81, 85] which can be tested in proton decay experiments in the 
near future. 


6 BAU and DM From Additional Leptoquarks 

As was seen in the previous section, having two light leptoquarks can be motivated by 
grand unification. If we now include an additional leptoquark in the model discussed 
above, we find additional freedoms for generating the baryon asymmetry. In fact, the 
gauge symmetries allow the operator 

Vp=q ijk S i S j S k H , (23) 

where i,j,k = 1,2. Although this operator breaks the discrete B — L discussed above, 
it may still give an acceptably small proton decay because the proton decay operators 
are induced at the two-loop level and the couplings between the leptoquark and the SM 
particles are extremely small. Another possible way to prevent proton decay would be to 
give the second leptoquark a different charge under B — L of 1/3. With this symmetry the 
above operator would still be allowed under the discreet B — L, but the second leptoquark 
would be forbidden from coupling to the SM particles alleviating proton decay. 

This operator can be relevant for generating the baryon asymmetry. Since the lepto¬ 
quarks are scalars, it is possible that during inflation one of the leptoquarks obtains a large 
VEV. After inflation, the VEV of the leptoquark will eventually return to zero. However, 
since the couplings can violate CP, the relaxation can generate angular momentum 
in the leptoquark fields [86]. This angular momentum violates baryon number and lepton 
number. The baryon number stored in the leptoquark angular momentum generates a 
baryon asymmetry in the quarks as the VEV decays. If conditions are right, this could be 
how the baryon asymmetry was generated. The lepton number stored in the VEV leads 
to the production of right-handed neutrinos when the VEV of the leptoquark decays. 
This could be a non-thermal means of producing the lightest right-handed neutrino. If 
the density is large enough, it could be DM. 

7 Conclusion 

There are several unexplained signals lurking at the LHC. One particularly interesting 
signal is the excess in both eejj and eujj found in the CMS first-generation leptoquark 
searches [8] . The signals found in these searches rely on the leptoquark having some decay 
paths which are hidden so the branching fractions that produce the eejj and eujj are 
sufficiently small. If we assume the leptoquark has SM charges (3, 2,1/6) in the basis 
SU(3) C <g) SU(2 ) l <8 ) U(l) y , it will have the interactions SdL and SQv R with generation 
indices suppressed. For this model, the simplest way of explaining the hidden decays is to 
assume that the leptoquark couples to all three generations. Since the current searches for 
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leptoquarks focus on leptoquarks that couple to each generation individually, the hidden 
decay modes could merely be decays to other generations. If this is the case, events with 
mixed-generation leptons like rejj should be seen at the upgraded LHC. 

On the other hand, if right-handed neutrinos are lighter than the leptoquark, they 
could decay to right-handed neutrinos. If the two heavier right-handed neutrinos are 
somewhat degenerate in mass with the leptoquark, these decay modes would be excluded 
because the jets produced in this event would be too soft. 

For each of these scenarios, the lightest right-handed neutrino is a viable DM candidate 
if its couplings to the leptoquarks are small enough that its lifetime is longer than the 
age of the Universe. If the coupling of the right-handed neutrinos to the leptoquarks are 
small enough, this model looks very similar to the ffiVISM [63] with similar production 
modes for the DM and baryon asymmetry. However, if these couplings are on the larger 
end of the allowed range and the reheat temperature of the Universe is small enough, 
the lightest right-handed neutrino DM can be a non-equilibrium thermal DM candidate 
[30, 78], 

With the leptoquark couplings small, the width of the leptoquark is also small. Thus, 
we expect a narrow peak in the invariant mass distribution constructed from a quark 
and a lepton final state. Since the statistics are low for this signal, we currently cannot 
conclude whether the excesses observed in the CMS experiment are actually peaked or 
not. Future LHC experiments with more data will be able to confirm or exclude this 
prediction. 

The quality of the gauge coupling unification in the SM is made more precise by the 
addition of a leptoquark. It is even better with two leptoquarks and an additional Higgs 
multiplet. Since proton decay in these SM extensions is governed by interaction with the 
heavy gauge bosons, the relative sizes of the different decay modes are determined by 
the CKM matrix elements. Since the mass of the heavy gauge bosons also tends to be 
relatively light, these models make fairly precise predictions for proton decay experiments. 

The additional leptoquarks also give other possibilities for producing the baryon asym¬ 
metry of the Universe. They allow for an Affleck-Dine [86] like production of the baryon 
asymmetry. 


Note Added: As we were about to post, we noticed Ref. [87] which also discusses an 
interpretation of the CMS excess [8] in terms of leptoquarks. However, this model is 
quite different from our model. They consider the leptoquark So instead of 5 T [ /2 which 
we consider. They also consider non-diagonal Yukawa couplings. The diagonality of onr 
couplings is advantageous because it does not have difficulty with rare meson decays. 
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